arXiv:1509.02856vl [astro-ph.HE] 9 Sep 2015 


ApJ, in press 

Preprint typeset using style emulateapj v. 5/2/11 


GAMMA-RAY FLARE ACTIVITY FROM PSR BI259-63 DURING 2014 PERIASTRON PASSAGE AND 

COMPARISON TO ITS 2010 PASSAGE 

G. A. CALIANDRO*'^, C. C. CHEUNG^ J. LP, J. D. SCARGLE^ D. F. TORRES"'-*, K. S. WOOd\ M. CHERNYAKOVA^’* 

ApJ, in press 

ABSTRACT 

PSR B1259-63/LS 2883 is a gamma-ray binary system containing a radio pulsar in a highly elliptical ~ 
3.4-year orbit around a Be star. In its 2010 periastron passage, multiwavelength emission from radio to TeV 
was observed, as well as an unexpected GeV flare measured by the Fermi Large Area Telescope (LAT). Here, 
we report the results of LAT monitoring of PSR B1259-63 during its most recent 2014 periastron passage. We 
compare the gamma-ray behavior in this periastron with the former in 2010 and And that PSR B1259-63 shows 
a recurrent GeV flare. The similarities and differences in the phenomenology of both periastron passages are 
discussed. 

Subject headings: gamma rays; stars — pulsars; individual (PSR B1259-63) — X-rays; binaries 


1. INTRODUCTION 

Gamma-ray binaries constitute a small subgroup of X- 
ray binaries hosting O/B stellar companions (Dubus 2013). 
They emit modulated radiation at the orbital period. In fact, 
the gamma-ray emission from these systems dominates their 
spectral energy distributions. Among all members in this 
group, PSR B1259-63 is the only one for which the compact 
object’s nature is known; a 47.76 ms radio pulsar orbiting a Be 
star (LS 2883) with a period of ^ 1236.7 days. It is located at 
a distance of 2.3 kpc (Johnston et al. 1992, 1994; Negueruela 
et al. 2011; Shannon et al. 2014). Because of the high incli¬ 
nation angle {io) between the pulsar orbital plane and the Be 
star circumstellar disk (in ^ 10° — 40°, Melatos et al. 1995), 
the circumstellar disk is crossed twice when the pulsar moves 
in one full highly elliptical (e ~ 0.87) orbit. During these 
crossings, multiwavelength emission results from the interac¬ 
tion between the relativistic pulsar wind and the stellar wind 
of the companion (Chernyakova et al. 2006,2009,2014; Aha- 
ronian et al. 2009; Abdo et al. 2011). 

Around the 2010 periastron passage of PSR B1259-63 
(2010 December 14), enhanced gamma-ray emission as well 
as a month-long GeV flare were observed (Abdo et al. 2011; 
Tam et al. 2011). During the GeV flare, PSR B1259- 
63 was characterized by an extremely high conversion effi¬ 
ciency of the pulsar spin-down power into gamma-ray emis¬ 
sion. The physical origin of the flare is still under debate (e.g., 
Khangulyan et al. 2012; Kong et al. 2012; Dubus & Cerutti 
2013; Chernyakova et al. 2014). After its 3.4-year orbit pe¬ 
riod, PSR B1259-63 approached periastron again (2014 May 
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4). We now report the observation of a flare following the 
2014 periastron at a time compatible with the extrapolated ex¬ 
pectation following the 2010 periastron. The initiation of the 
flare was reported in Astronomer’s Telegram^ by Malyshev 
et al. (2014); Tam et al. (2014); Wood et al. (2014). Detailed 
results of our monitoring campaign in GeV and X-rays, ana¬ 
lyzing the Fermi Large Area Telescope (LAT; Atwood et al. 
2009) data and SwiftlXKY data, respectively, are given. We 
also provide a comparison with a reanalysis of the 2010 peri¬ 
astron data using the latest LAT instrument response functions 
(IRE) and diffuse models H 

2. OBSERVATIONS AND DATA ANALYSIS 

The analysis of Fermi-LAT data was performed using the 
Fermi Science ToolJ^ 09-34-01 release. For the 2010 and 
2014 periastron passages the analysis was carried out with 
Pass 7 reprocessed data belonging to the SOURCE event 

clas£3- 

All gamma-ray photons within an energy range of 0.1-100 
GeV and within a circular region of interest (ROI) of 10° ra¬ 
dius centered on PSR B1259-63 were used for this analysis. 
To reject contaminating gamma rays from the Earth’s limb, 
we selected events with zenith angle < 100°. Fermi-LAT had 
triggered Target of Opportunity (ToO) observations in 2014 
starting 27 days after the periastron, and lasting 19 days. Dur¬ 
ing the ToO, the exposure of PSR B1259-63 was increased 
^ a factor of ^2 compared to normal observations (Figure 
[T] panel b). In contrast, during the 2010 periastron passage 
a modified sky survey had commenced 13 days after perias¬ 
tron. In this mode, which lasted 10 days, the southern hemi¬ 
sphere receives 30% extra exposure (Figure [T] panel a). The 
gamma-ray flux and spectral results of PSR B1259-63 pre¬ 
sented in this work were calculated by performing a binned 
maximum likelihood fit using the Science Tool gtlike. The 
spectral-spatial model constructed to perform the likelihood 

'^http: //www.astronomerstelegram.org/ 

***111 the process of preparing this paper, a similar study by Tam et al. 
(2015) was presented. While their main results are consistent with ours, we 
have included a more detailed comparison between 2010 and 2014 flares as 
well as pulsar timing analysis. 

"http://fermi.gsfc.nasa.gov/ssc/ 

*^http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 
documentation/Pass7REP_usage.html 
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Fig. 1.— Fermi-LAT exposure of PSR B1259-63 during the 2010 (panel 
a) and 2014 (panel b) periastron passages. 

analysis includes Galactic and isotropic diffuse emission com¬ 
ponents as well as known gamma-ray sources within 15° of 
PSR B1259-63 based on a catalog internal to the collabora¬ 
tion (now released as the 3FGL catalog; Acero et al. 2015). 
The spectral parameters were fixed to the catalog values, ex¬ 
cept for the sources within 3° of PSR B1259-6S For these 
latter sources, the flux normalization was left free. We veri¬ 
fied that changing 3° to 5° produced consistent results. PSR 
B1259-63 itself was modeled as a single power-law with all 
spectral parameters allowed to vary. The contribution of the 
Galactic and isotropic diffuse emissions within the analyzed 
ROI was estimated performing a preliminary maximum like¬ 
lihood fit. For the 2010 periastron event, this fit included data 
from 5 months prior to 1 year after periastron. For 2014, 
this fit included data from 5 months prior to periastron. The 
resulting scale factors of the isotropic and Galactic diffuse 
templates were kept fixed in the analysis that follows. The 
2010 and 2014 periastron times used in this paper are MJD 
55544.693781 (2010 December 14 16:39:02.678) and MJD 
56781.418307 (2014 May 4 10:02:21.725), which were de¬ 
rived from the most recent measurement of the binary sys¬ 
tem’s orbital parameters (Shannon et al. 2014). 

SwiftfKKV also monitored PSR B1259-63 during the 2014 
periastron passage in Photon Counting (PC) and Window 
Timing (WT) modeS We selected PC data with event grades 
0-12 and WT data with event grades 0-2 (Burrows et al. 
2005). Source events were accumulated within a circular re¬ 
gion with a radius of 30 pixels (1 pixel = 2.36 arcsec). Back¬ 
ground events were accumulated within a circular, source-free 
region with a radius of 60 pixels. Exposure maps were gen¬ 
erated with the task XRTEXPOMAP. Ancillary response flies 
were generated with XRTMKARE, which accounts for differ¬ 
ent extraction regions, vignetting, and point spread function 
corrections. We analyzed the SwiftlXKY 0.1-10 keV data by 
HEAsoft version 6.14,0 and spectral fitting was performed 
using XSPECV. 12.8.1. 

3. LIGHT CURVE 
3.1. GeVflares 

*^http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 
scitools/source_models.html 

^^https : / /www.swift.psu. edu/xrt / software .html 
*^http:/ /heasarc . nasa. gov/lheasoft/ 


During PSR B1259-63’s 2010 periastron passage, bright 
GeV emission was observed ~ 30-80 days after periastron 
(Abdo et al. 2011). In the 2014 periastron passage, we moni¬ 
tored the gamma-ray emission from PSR B1259-63 on daily 
and weekly timescales starting 60 days before periastron. Eor 
the time periods when PSR B1259-63 was not detected with 
the EAT (Test Statistic (TS0 value < 9), we placed a 95% 
confidence upper limit on the photon flux above 100 MeV, 
as evaluated with Helene’s method (Helene 1983). A rapid 
GeV flux and TS value increase was observed beginning 2014 
June 6 (31 days after periastron, Eigure|2] red points). The 
GeV flare is apparent to the eye when plotted in a weekly 
timescale (Eigure lUt). The detailed daily light curve of the 
2014 flare is shown in Eigure |2ji, |2) &|2]. A large variability 
is observed in the 0.1-100 GeV band, with multiple peaks in 
a daily timescale. To compare the 2010 flare and 2014 flare, 
we reanalyzed the data from the 2010 periastron passage. The 
weekly and daily light curves showed consistent results with 
those reported in Abdo et al. (2011), in which the analysis 
was carried out with Pass 6 data and P6_V3 IRlFI. 

In spite of obvious similarities, the two flares show differ¬ 
ent time evolution. The light curves of the 2010 (blue) and 
2014 (red) periastron passage are plotted together in Eigure 
|2l It is apparent to the eye that in the weekly light curve, the 
2010 flare has a higher peak flux and a rapid decreasing evolu¬ 
tion (Eigure|2^). The 2014 flare has a lower peak flux, which 
then persists rather than falling rapidly. The daily light curves 
show these same trends. The 2010 flare consists of a large 
outburst in the beginning followed by several weaker flares. 
The 2014 flare is composed of several sub-flares with similar 
flux levels. 

We calculated the cross-correlation using crossco^ in 
HEASOET with the daily flux of the 2010 and 2014 flares 
(30 days to 80 days after periastron), which reflects the differ¬ 
ence of general trend. We found that the 2014 flare is 2.4±0.6 
days delayed from the 2010 flare. The error is estimated with 
10000 simulations by sampling the daily light curve fluxes 
within the errors and calculating the cross-correlation for each 
one. 

In order to better visualize similarities and differences be¬ 
tween the 2014 and 2010 flare profiles, their smoothed light 
curves are plotted in Eigure [3 The smoothed light curves 
were produced using a sliding window technique. We chose 
time windows of 3 days, whose starting times lag the previous 
one by 3 hours. Eor each time window, a binned likelihood 
analysis was performed. The spectral index of PSR B1259- 
63 was allowed to vary between the values 2.0 and 3.5. With 
the sliding window technique the smoothing is due to the cor¬ 
relation of the points at a time distance smaller than the time 
window. 

Eigure [^ & [3i show the flux and the TS smoothed light 
curves of the 2010 and 2014 periastron passages. In Eigure]^ 
the shaded areas around the solid lines correspond to the sta- 

*®The Test Statistic is defined as TS=—2 ln(Lmaa:,o/^maa;,l), where 
Lmaxfi is the maximum likelihood value for a model without an additional 
source (the “null hypothesis”) and Lmax,l is the maximum likelihood value 
for a model with the additional source at a specified location. A larger TS 
indicates that the null hypothesis is inconect (i.e., a source really is present). 
As a basic rule of thumb, the square root of the TS is approximately equal to 
the detection significance for a given source. 

*'^http; // fermi . gsfc.nasa . gov/ssc/data/analysis/ 
documentation/Cicerone/Cicerone_LAT_IRFs/ 
IRF_overview.html 

**https://heasarc.gsfc.nasa.gov/xanadu/xronos/ 
help/crosscor.html 
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Fig. 2.— Weekly (left panels) and daily (right panels) gamma-ray flux (a & b, g & h), TS value (c & d, i & j), and photon index (e & f, k & 1) of PSR B1259-63 
during the 2010 (blue) and 2014 (red) periastron passages. The dashed black line in left panels is a fiducial zero indicating the time of periastron. The shaded 
yellow and green regions indicated in panels g& h are the defined Peak and Tail periods. 


tistical error of the flux calculation. When the TS value is low 
(TS < 9), the smoothed light curve points are represented as 
null flux with shaded areas showing the upper limits of 95% 
confidence. In the TS smoothed light curves, the onset of 
the flaring activity in 2010 and 2014 periastron start approx¬ 
imately at the same time, ~ 31 days after periastron. After 
the onset, the 2010 periastron smoothed flux steeply increases 
up to its maximum at ^ 36 days. However, in 2014 the in¬ 
crease is delayed by forming a short plateau until ^ 34 days, 
after which the flux rises up to the peak (at ^ 38 days). The 
main peak of the 2014 profile is delayed by ~ 2 days with 
respect to the 2010 flare, which is consistent with the result of 
flux cross-correlation. A further significant difference is that 
in 2010 the main peak is followed by a marked second peak 
(at r\j 43 days), which in the 2014 light curve is substituted by 
a smooth valley. Finally, the 2014 flaring activity ends with a 
possible peak at ~ 68 days, which is consistent with the last 
high flux bin of the weekly light curve and high upper limits 
of the daily light curve (see Figure[2]i. However, the daily ex¬ 
posure time of this possible peak is ^ 5 times less than the 
main peak, which leads to low TS values and makes this peak 
insignificant, preventing further analysis. 

We searched for gamma-ray pulsations from PSR B1259- 
63 following the same weighted H-test procedures described 
in Abdo et al. (2011) and Hou et al. (2014), using a radio 
ephemeris internal to the collaboration measured at Parkes 
Observatory. All LAT observations available excluding pe¬ 
riastron passages (from 60 days before periastron until 120 
days after periastron) were taken into consideration. Neither 


a statistically significant detection of a pulsation, nor a detec¬ 
tion of the source was found. We also searched for gamma-ray 
pulsation around the 2009 and 2012 apastron passages (from 
2,4 and 6 month intervals centered on apastron). Neither pul¬ 
sation nor the source itself were significantly detected in these 
periods. Gamma-ray pulsations in the timeframe around the 
2010 and 2014 periastron passages (from 60 days before peri¬ 
astron until 120 days after periastron) were also searched for, 
but without a detection. 

3.2. Multiwavelength view of both periastra 

PSR B1259-63 was monitored by 5w!/f/XRT during its 
2014 periastron. No significant flux increase on timescales 
shorter than ^4 days was observed in soft X-rays during the 
GeV flare (Chernyakova et al. 2015). The light curve follows 
the general trend of previous periastron measurements (Fig¬ 
ure®?). X-ray measurements of PSR B1259-63 during 2014 
and previous periastra, and radio measurements during 2010 
and 1997 periastra all show a double peaked structure (Fig¬ 
ure lU? & |4};), which is expected when the pulsar crosses the 
companion’s circumstellar disk twice each orbit (Cominsky et 
al. 1994; Chernyakova et al. 2006, 2009; Abdo et al. 2011). 
In Figure |4]the position of the Be circumstellar disk is shown 
with shaded areas as proposed in Chernyakova et al. (2006) 

In the 2010 and 2014 first crossings of the Be circumstellar 
disk, the gamma-ray flux measurements are upper limits be¬ 
cause of the low TS value. Although the 2010 and 2014 flares 
are within the second crossings of the Be circumstellar disk, 
they are not in line with X-ray and radio activities (Figure|4^). 
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Fig. 3.— Flux (a) and TS (b) smoothed light curves of the 2010 and 2014 
flaring activities after periastron (blue and red, respectively) produced with 
a sliding window. The shaded areas demonstrate the statistical eiTor of the 
flux. We note that we removed from the plot a few data points where the 
maximum-likelihood fit did not provide reliable error bars. The smoothed 
light curve points are represented as null flux when the TS value is below 9, 
with only shaded areas showing the upper limits at 95% confidence. See text 
for more details. 


TABLE 1 

Power-law spectrum parameters of PSR B 1259-63 during 2014 
AND 2010 FLARES 




Photon index 

r 

FluxOlOOMeV) 
10“’^ ph cm-2 s-l 


Average 

2.93 ± 0.07 

9.8 ±0.7 

2014 flai-e 

Peak interval 

2.99 ±0.10 

II.3 ± I.O 


Tail interval 

2.87 ± 0.09 

8.6 ±0.9 


Average 

2.90 ± 0.07 

II.8±0.9 

2010 flare 

Peak interval 

2.98 ±0.10 

18.8 ± 1.5 


Tail interval 

2.83 ± 0.09 

8.8 ± 1.1 


3.3. Spectral analysis 

On the basis of the smoothed light curves we defined two 
time intervals of the flare profiles. The first interval includes 
the initial peak (labeled “Peak” interval). All the rest of the 
flare constitutes the second interval (labeled “Tail” interval). 
The transitions from Peak to Tail are set when TS values in 
Figure [3 reach the lowest value soon after the peak. For the 
2010 flare, the transition from Peak to Tail occurred 40 days 
after the periastron while for the 2014 flare this transition oc¬ 
curred 42 days after periastron (Figure |2^ & h). The begin¬ 
nings of both flares occur at 31 days, which is the start of the 
TS value & flux rising in Figured The ends of both flares oc¬ 
cur at 71 days, which is the time when the TS values in Figure 
^finally drop below 10. 

The average spectra of the whole flares, as well as the spec¬ 
tra of the Peak and Tail intervals were evaluated performing 
a binned likelihood analysis. The results are shown in Table 
1 and Figure |5]. The average spectra during the 2010 and 
2014 GeV flare periods are modeled by a single power law 
and a power law with an exponential cutoff. The spectra are 
shown in Figure We compare the two models utilizing the 
likelihood ratio test (Mattox et al. 1996). The A TS between 
the two models is less than 9, which indicates that the signifi¬ 


cance of the spectral cutoffs for 2010 and 2014 flares are less 
than 3(7. Thus, a power law with an exponential cutoff does 
not significantly improve the fitting over the simpler power- 
law model. The fitted parameters of a single power law model 
are shown in Table 1. The two flares show comparable aver¬ 
age flux levels and spectral shape (Table 1; Figure [Sjt &|5]:). 
No spectral variability is seen above 3a during both 2010 and 
2014 flares. The flux of the Peak interval in the 2010 flare is 
significantly higher than in the 2014 flare (Figure l^t), while 
the Tail intervals have similar fluxes (Figure|2;). In both flares 
the spectral slope is steeper in the Peak interval than in the 
Tail, but the variations are not significant. 

4. DISCUSSION 

We have observed the GeV flare during the 2014 perias¬ 
tron passage of PSR B1259-63 with Femu'-LAT These are 
our main results; 

• We confirm a recurrent flaring behavior in gamma rays, 
a phenomenon that is associated with the periastron 
passages. 

• The 2014 GeV flare exhibits a similar average flux level 
and spectral shape with the 2010 flare. 

• The two GeV flares showed different flux evolution. 

• The 2014 GeV flare peak is about 2 days delayed from 
the 2010 flare. 

• No GeV pulsations from PSR B1259-63 have been de¬ 
tected in any part of its orbit. 

We have shown that the 2010 and 2014 periastron passages 
of PSR B1259-63 produces some striking similarities in their 
gamma-ray emission. At the same time, there are also certain 
differences that may be due to inhomogeneities in the shape, 
density, or extent in the circumstellar disk of the Be star. The 
two flares were characterized by a high efficiency in the con¬ 
version of pulsar spin-down power into gamma-ray emissions. 
The 2014 periastron highest day-average flux was 7.1 x 10“^° 
erg cm“^ s“^ with a photon index of ^^3.1. Though it is 
lower than the peak day-average flux of the 2010 periastron, 
9.6 X 10“^° erg cm“^ s“^, it corresponds to an isotropic 
gamma-ray luminosity of 4.5 x 10^® erg s“^, equaling ^ 50% 
of the pulsar spin-down luminosity, 8.3 x 10^^ erg s“^ (John¬ 
ston et al. 1992), without considering possible beaming ef¬ 
fects. Gamma-ray conversion efficiency of most gamma-ray 
emitting pulsars is below 30% (Abdo et al. 2013). Though we 
have now observed two similar GeV flares from PSR B1259- 
63, the origin of these events is still unclear. Khangulyan et 
al. (2012) account for the GeV flare by the upscattering of the 
circumstellar disk IR photons by the unshocked pulsar wind 
electrons, which will increase after the pulsar exits the sec¬ 
ond crossing of the disk. Kong et al. (2012) explain the GeV 
flare as the Doppler boosted synchrotron emission in the bow- 
shock tail, which could accommodate the high conversion ef¬ 
ficiency of pulsar spin-down power. Dubus & Cerutti (2013) 
adopted upscattering of X-ray photons from the pulsar wind 
nebula, which would produce correlated X-ray and gamma- 
ray emission. Such correlation was not observed, which could 
be due to Doppler boosting or a realistic scattering geometry 
of the system. The unresolved nature and competing models 
of GeV flares will motivate future observations dedicated to 
PSR B1259-63 and hopefully improve our understanding of 
its unique phenomenology. 
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Fig. 4.— Light curves of PSR B1259-63 around periastron. (a) weekly Fermt-LAT gamma-ray flux of 2010 (blue) and 2014 (red) periastron passage in energy 
units, (b) X-ray fluxes from 2014 Swift (red) and previous periastron passages (blue, XMM-Newton, Chandra, Swift and Suzaku observations during 2004, 
2007 and 2010 periastra, Chernyakova et al. 2009. The errors are smaller than the symbols), (c) radio (2.4 GHz) flux densities for the 2010 (blue) and 1997 
periastron passages (green, Abdo et al. 2011). X axis is in true anomaly. The dashed black line shows the time of periastron. Shaded area corresponds to the Be 
circumstellar disk position proposed in Chernyakova et al. (2006). 
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Fig. 5.— F(erm/-LAT flare spectra of 2010 (blue) and 2014 (red) periastra. 
(a) Spectra calculated over the whole duration of the flares, (b) Spectra of 
the Peak intervals of 2010 and 2014 flai'e. (c) Spectra of the Tail intervals of 
2010 and 2014 flare. 
































































